Computational ghost imaging versus imaging laser radar for 3D imaging 
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Ghost imaging has been receiving increasing interest for possible use as a remote-sensing system. 
There has been little comparison, however, between ghost imaging and the imaging laser radars 
with which it would be competing. Toward that end, this paper presents a performance comparison 
between a pulsed, computational ghost imager and a pulsed, floodlight-illumination imaging laser 
radar. Both are considered for range-resolving (3D) imaging of a collection of rough-surfaced objects 
at standoff ranges in the presence of atmospheric turbulence. Their spatial resolutions and signal- 
to-noise ratios are evaluated as functions of the system parameters, and these results are used to 
assess each system's performance trade-offs. Scenarios in which a reflective ghost-imaging system 
has advantages over a laser radar are identified. 

PACS numbers: 42.30.Va, 42.68.Bz, 42.68.Wt 



I. INTRODUCTION 

Ghost imaging is an active-imaging technique that 
uses time- varying structured illumination to image a tar- 
get without spatially-resolving measurements of the light 
beam that interacts with the target. Traditionally, a 
beam splitter is used to create two perfectly-correlated 
beams, the signal and reference, such that the signal in- 
teracts with the target and is then measured by a single- 
pixel bucket detector, while the reference is directly mea- 
sured by a spatially- resolving detector As the illu- 
mination pattern is varied, the two measurements are 
correlated until the spatial structure of the target is de- 
termined. Neither measurement alone is sufficient to pro- 
duce the image; it is their cross-correlation that holds the 
desired target information, i.e., the ghost image. 

The first ghost-imaging experiment relied on entan- 
gled photon-pairs — obtained from spontaneous paramet- 
ric downconversion — for its signal and reference fields, so 
it was believed that the ghost image was a uniquely quan- 
tum feature 0], namely nonlocal two-photon interfer- 
ence. Later, ghost imaging was performed with classical 
pseudothermal light 0,13, and a controversy arose as to 
whether these experiments could be explained by the in- 
tensity correlation between classical signal and reference 
fields or only by nonlocal two-photon interference. Sub- 
sequent Gaussian-state analysis provided a unified treat- 
ment of downconverter and pseudothermal ghost imag- 
ing, showing that the stronger-than-classical correlation 
of entangled photons yielded better contrast and near- 
field resolution [H, H, [|| , and that the quantum and semi- 
classical treatments of the pseudothermal imager gave 
quantitatively-identical performance predictions. More 
recently, the intensity cross-correlation versus nonlocal 
two-photon interference controversy for understanding 
pseudothermal ghost imaging has been ended through 
analysis demonstrating that these two explanations can 
co-exist Q. 

Once pseudothermal ghost imaging is considered in 
the framework of structured-illumination imaging, it be- 
comes possible to dispense with a physically-realized ref- 



erence field. In particular, deterministic modulation of a 
laser beam with a spatial light modulator (SLM) can pro- 
vide the signal field used for target interrogation, while 
the on-target intensity pattern needed for the reference 
field can then be calculated via diffraction theory [8J. 
The extension of ghost imaging to this computational 
framework has opened the door for a variety of appli- 
cations, including demonstration of ghost imaging with 
phase-sensitive classical light Q, and image reconstruc- 
tion via compressed sensing, instead of correlation [Toj . 
Pseudothermal ghost imaging has been experimentally 
done in reflection both in laboratory 111 and remote- 
sensing scenarios Moreover, recent work analyzing 
reflective ghost imaging indicates that computational re- 
flective ghost imaging is feasible for remote sensing [l3j |. 
and computed reference beams can be generated for all 
target ranges of interest, so that computational ghost 
imaging has unlimited depth of focus, unlike its pseu- 
dothermal counterpart. So far, however, there has only 
been a cursory performance comparison between compu- 
tational ghost i mag in g and an imaging laser radar for 
this application [HI. fl5|. 

In this paper we extend the analysis in jl3| of a 
continuous- wave reflective ghost imager to a pulsed, com- 
putational ghost-imaging system that is capable of per- 
forming three-dimensional imaging, and we compare its 
behavior to that of a pulsed, floodlight-illumination, 
imaging laser radar. We present results for their spatial 
resolutions and signal-to-noise ratios (SNRs) when imag- 
ing rough-surfaced targets that produce fully-developed 
laser speckle and the propagation to and from the tar- 
gets is through atmospheric turbulence. We also in- 
vestigate the trade-off between spatial resolution and 
SNR as a function of detector and entrance-pupil sizes. 
As reflective ghost-imaging systems are only subject to 
spatial-resolution loss from turbulence in the source-to- 
target path [l3|, [HI, whereas a floodlight-illumination 
laser radar's spatial resolution is only degraded by turbu- 
lence in the target-to-receiver path, we consider arbitrary 
turbulence distributions on all optical paths. 
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FIG. 1: (Color online) (a) Setup for 3D computational ghost 
imaging in reflection. Pulsed laser-light undergoes spatial 
light modulation, propagation through atmospheric turbu- 
lence to and from an extended range-spread target, and shot- 
noise limited bucket detection, producing an output, if,(n, z), 
from the range- z return associated with the nth transmitted 
pulse. Diffraction theory is used to calculate I(p',n,z), the 
vacuum-propagation, target-region intensity pattern at trans- 
verse coordinate p' and range z associated with the nth trans- 
mitted spatial pattern, (b) Setup for 3D-imaging laser radar. 
Pulsed laser-light illuminates the target region, in a floodlight 
manner, through a turbulent optical path. The return light, 
collected after propagation back to the radar through atmo- 
spheric turbulence, is detected by a shot-noise limited CCD 
array in an image plane, resulting in an output, i p (p p ,n, z), 
for the pixel at transverse coordinate p' p and range z from the 
nth transmitted pulse. 



II. SOURCE, PROPAGATION, AND TARGETS 

Figures 1 (a) and (b) show the setups that we shall con- 
sider for computational ghost imaging and laser radar, 
respectively. In this section we highlight the features 
that are common and different in these systems. Their 
performance characteristics will be treated in Sees. Ill 
(ghost imaging) and IV (laser radar). For both imagers, 
a pulse-train waveform is emitted by the source and prop- 
agates through atmospheric turbulence to an extended, 
range-spread, rough-surfaced target. The light reflected 
from the target propagates back through atmospheric 
turbulence to the receiver, where it is photodetected. For 
the ghost imager, each pulse has a different spatial pat- 
tern, whereas the laser radar uses a constant, floodlight- 
illumination pattern for all its pulses. The ghost im- 
ager employs a single-pixel bucket detector, and forms 
its image by cross-correlation with a computed intensity- 
pattern reference. The laser radar, however, employs a 



focusing lens and a CCD array to form an image of the 
target region. To make a fair performance comparison 
between these two systems, we will require that their 
sources radiate the same average photon number per 
pulse, Nt, towards the target region. In the subsections 
that follow, we provide source, propagation, and target 
details that will serve as the foundation for the perfor- 
mance analyses to come in Sees. Ill and IV. 



A. Source Characterization 

We will take the SLM output in Fig. 1(a) and the laser 
output in Fig. 1(b) to be quasimonochromatic, classi- 
cal scalar waves with center frequency luq (wavelength 
Ao = 2ttc/ujo, wave number fco = ujq/c) and complex en- 
velope E s (p, t) as a function of transverse coordinate, p, 
in the transmitter's exit pupil and time, t. We normal- 
ize E s (p,t) so that \E s (p,t)\ 2 is the photon flux-density 
emitted from p at time t. The pulse-train waveforms for 
the ghost imager and the laser radar that we will employ 
are 



N-l 



E s (p,t) = ^^2^,n)p(t-nT B ) t 



(1) 



71=0 



where: N is the number of pulses to be used in forming 
the image [l7j]; n) is the normalized (/ dp |£(p, n)\ 2 = 
1) spatial mode of the nth pulse; p(t) is a normalized 
(Jdt\p(t)\ 2 = 1) pulse shape that is time-limited to 
\t\ < Tp/2; and T s > T p is the pulse- repetition inter- 
val. The pulse duration T p , which could be picoseconds 
to nanoseconds long, will be taken to be much shorter 
than the coherence time, T c , of the turbulence, which is 
typically milliseconds long. The pulse-repetition interval, 
T a , will be long enough to preclude second-time-around 
echoes from the farthest range of interest to masquer- 
ade as returns from a closer range [l8j]. Also, T s will be 
allowed to exceed T c . 

The principal difference between the fields transmit- 
ted by the ghost imager and the laser radar lies in their 
normalized spatial modes, £(p, n). For the ghost imager 
these will be modeled as a collection of independent, iden- 
tically distributed, zero-mean Gaussian random fields — 
indexed by the pulse number n — that are completely 
characterized by the following Gaussian Schell-model cor- 
relation functions fijl 



(t(Pi,ni)£(p2,n 2 )) = 0, 



(2) 



-(\Pi\ 2 + \P2\ 2 )/a 2 



(3) 

Here, ao is the source's intensity radius and po "C ao 
is its spatial coherence length. The Gaussian Schell- 
model is an analytical convenience that captures the 
essential physics of radiation from a partially-coherent 
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source. The reader is cautioned, however, that the Gaus- 
sian statistics we have assumed will not be valid until 
sufficient propagation away from the SLM has occurred 
that the Central Limit Theorem can be applied to the su- 
perposition of the many phase-modulated field elements 
from that modulator's pixels. Inasmuch as we will be 
taking propagation from a single SLM pixel to be in its 
far-field regime, the Gaussian assumption will indeed be 
applicable for the ghost imager's target illumination. 

The normalized spatial mode for the laser radar will 
be deterministic, identical for all pulses, and given by 
the collimated Gaussian beam 



£(p.«) = 



nwi 



o-\p?h 



where w is its beam waist. For a fair comparison with 
the ghost imager, we shall set the laser radar's beam 
waist w equal to the ghost imager's p a . By doing so, 
the average far-field intensity pattern produced by the 
latter's source spatial-mode will exactly match the far- 
field intensity pattern produced by that of the former. 



B. Propagation through Turbulence 

Propagation of E s (p,t) from the ghost imager or the 
laser radar to range z in the target region — and propa- 
gation of a reflected field E r (p',t) from range z back to 
those sensors — will be through atmospheric turbulence, 
whose behavior we will characterize via the extended 
Huygens-Fresnel principle (20l |. Labeling the path from 
the sources to the target region with S, and the path 
leading to the detectors with D, we have that 



fluctuations imposed on the field arriving z m downrange 
at time t and transverse coordinate p' from a point 
source at transverse coordinate p; there is a correspond- 
ing physical interpretation for the real and imaginary 
parts, XD(p,p\t) and (f) D (p, p' ,t), of i/j d (p, p', *)• 

The range-spread targets of interest will be assumed 
to lie at ranges between z m i n and z max = z m ; n + Az from 
the two imagers under consideration, where -z m i n satisfies 
the far-field propagation conditions, kQaopo/2z m - m <C 1 
and fco?A>o/2z m i n *C 1, respectively, for the ghost im- 
ager's correlation functions, and laser radar's transmit- 
ter beams. As will be explained below — in conjunction 
with our spatial-resolution analyses for the ghost imager 
and the laser radar — we will not assume that far-field 
propagation conditions automatically apply for the ghost 
imager's source diameter or the laser-radar receiver's en- 
trance pupil diameter. 

We will assume that the S and D paths are suffi- 
ciently separated — by virtue of our imagers having differ- 
ent transmitter exit optics and receiver entrance optics — 
that ips(p',p,t) and iPd(p, P! ,t) will be statistically 
independent. It then turns out that turbulence will 
only enter in our spatial-resolution calculations through 
the mutual coherence functions of cxp[V>s(p', p, t)] and 
exp[ipD{p,p\t)]. Given our assumption of Kolmogorov 
spectrum turbulence, the mutual coherence functions we 
will need are as follows, 



/gV's (p'i -Pi .* ) e 1>s (Pa ,P2 ,t ) \ 



-£ ) s(pi-P 2 iPl-P2)/2 



(9) 



and 



l e i'n(p-i->p'i> t ) e ^D{p2,p' 2 ,t)\ _ e -D D { Pl -p 2: p' 1 -p' 2 )/2 ^ /jq\ 

where 



E z (p', t)= dp E s (p, t - z/c)h s z (p', p, t), 



and 



where 



and 



K(P,p\t) 



^ O gjfco(z+|p-p' | 2 /2z) e ipD(p,P ,t) 



i2nz 



(5) 



E d (p,t)= JdpE r (p',t-z/c)h°(p,p',t), (6) 



kne iko (z+\p'-p\ 2 /2z) e V's (p' ,p,t) 

h S z (p\ P, t) = -5 (7) 

ilnz 



(8) 



specify the atmospheric propagation kernels at time t 
in terms of the Fresnel-diffraction Green's function and 
the prevailing atmospheric fluctuations, tps(p',p, t) and 
^>£)(p, p' \t), arising from Kolmogorov-spectrum turbu- 
lence distributed, in general, non-uniformly along two 
paths. The real and imaginary parts, xs(p' P,~t) arL d 
4>s(p' , P, t), of ips(p', p, t) are the logamplitudc and phase 



D s (p',p) = 2Mk 2 z / dsCl s (sz)\p's + p(l-s)\ 3 , 
Jo 

(11) 

and 



D D (p,p') = 2.91k 2 z / dsCl D (sz)\p(l-s)+p's\* , 
Jo 

(12) 

with { C* S (C) : < C < z } and { C* D (C) : < C < z } 
being the turbulence-strength profiles on the 5* and D 
paths from the imagers' location ( = to the target 
region at range z. In order to obtain closed-form results, 
however, we shall replace these 5/3-law mutual coherence 
functions with their square-law approximations, i.e., we 
will use Eqs. © and (JTUJ) with 



Ds(p\p) 
and 

D D {p,p') 



p'\ 2 W' s + p 1 - p(8/3 -W' s - W s ) + \p\ 2 W s 



Pi 



(13) 



\p\ 2 W D + p p'(8/3 -W D - W' d ) + \p'\ 2 W' D 



Pi 



(14) 
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In these equations, {W m , Wj^}, for m = S, D, are path- 
weighting terms, given by 



w 



ds(l - s) C„ m (sz) 



(sz) 



(15) 



(16) 



in terms of the normalized turbulence-strength profiles, 



^n,m(C) 



Cn,m(C) 



Io dsC n, m ( sz ) 



for m = S, 



(17) 



and {ps, Pd} are the spherical- wave turbulence coherence 
lengths for the S and D paths for uniform turbulence 
distributions with the same integrated strength as the 
actual distribution, 



Pm = (l.09h$zj dsC 



2 



(sz) 



-3/5 



for m = S, D. 



(18) 

In our treatment of ghost imaging we will assume that 
Pa <C ps/\/Ws, as will typically be the case, but we 
will allow ps/y/Ws to be larger or smaller than ao- For 
the laser radar, we have wq = po, so we can assume 
w o "C Ps/y/Ws, but we will allow Pd/VWd to be larger 
or smaller than ri, which is the radar-receiver's entrance 
pupil radius. 



transverse dimensions, in excess of our sensors' spatial- 
resolution capabilities. Then, the field returned to either 
the ghost imager or the laser radar is given by 



K 



E d { P ,t) = ^r 



k=l 



J dp' T Zk (p')E Zk (p' ,t 



z k /c)h°(p,p',t), 
(19) 

where E Zk (p 1 , t), the target illumination at transverse co- 
ordinate p' , range z k , and time t, is obtained from Eq. ([5]) 
with z = Zk, and T Zk (p') is the (random) field- reflection 
coefficient for the target at range Note that Eq. (TT9")) 
assumes that the reflector at a particular range does not 
occlude those that are farther away from the sensors. 

To complete our target-return model, we only need to 
supply statistics for the {T Zk (p')}. Following laser radar 
theory [221 ]. these will be taken to be statistically inde- 
pendent, zero- mean, complex- valued Gaussian random 
processes characterized by the autocorrelation function 



(20) 



where T Zk (p 1 ) is the average intensity-reflection coeffi- 
cient for the reflector at range zu- Strictly speaking, the 
Gaussian statistics cannot apply very close to the target, 
but, because of the quasi-Lambertian nature of rough- 
surface reflections, the Central Limit Theorem will ap- 
ply to Ed(p,t), so our assuming that the {T Zk (p')} have 
Gaussian statistics is indeed warranted. 



C. Target Reflection 

Spatially-coherent, quasimonochromatic light reflected 
by an opaque object whose surface is rough on the scale 
of the illumination's wavelength — as most real- world sur- 
faces are — yields speckles in the resulting far-field inten- 
sity profile [2l[ . The speckles arc due to the superposition 
of randomly phase-shifted reflections from surface facets 
possessing microscopic wavelength-scale height variations 
with correlation lengths on the order of a wavelength. 
Averaged over the speckle behavior, the rough surface 
is a quasi-Lambertian reflector that sends light back into 
the hemisphere with some average intensity-reflection co- 
efficient. For both our computational ghost imager and 
our laser radar, it is this intensity-reflection coefficient, 
as a function of transverse coordinate and range, that 
is to be imaged. In both cases, the speckle fluctuations 
present a significant impairment to high-quality image 
formation. Hence we need to incorporate a statistical 
model for speckle — relating the field illuminating an ex- 
tended range-spread target to the field returned to our 
sensors — into our analysis. 

As suggested in Fig. 1, we will presume the target to 
comprise a set of K quasi-planar, rough-surfaced objects 
that are located at discrete ranges, { Zk : 1 < k < K}, 
obeying z min < z\ < z 2 < ■ ■ ■ < z K < z max . Moreover, 
these objects will be taken to have spatial extents, in 



III. COMPUTATIONAL GHOST IMAGER 

A. Image Formation 

The computational ghost image for range z is formed 
by cross-correlating the computed intensity patterns, 
{I(p',n,z) : < n < N — 1}, that would be pro- 
duced from vacuum propagation of the TV transmitted 
spatial patterns {£(p, n) : < n < N — 1}, with the 
corresponding sequence of bucket-detector outputs that 
have been matched filtered for range z, which we denote 
{ib(n, z) :0 <n < N — 1}. The {I(p' ,n, z)} are given 
by 



I(p',n,z) 



k e 



ik \p -p\ 2 /2z 



i2nz 



(21) 



For convenience we have omitted leading constants that 
would appear if we used a physical reference; this scal- 
ing does not affect our analysis. The matched-filtered 
photocurrents {ib(n, z)} satisfy 



i b (n,z) 



J dr g(nT s + 2z/c- 



r) 



qr, dpA 2 b (p)\E d {p,T)\ 2 + Aift(T) 



,(22) 
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where: g(t) = \p(— t)\ 2 is the matched filter (causality 
ignored) for the transmitted pulse's intensity; q is the 
electron charge; rj is the bucket detector's quantum effi- 
ciency; Ab(p) = e~' p ' l T ' b is the field-transmission pupil 
function for the bucket detector's photosensitive region 
[23| ; and Aib(t) is the photocurrent shot noise, which we 
have assumed to be entirely due to the target return. (For 
a treatment of ghost imaging that includes background 
light and its associated shot noise see [H[.) 
The TV-pulse ghost image for range z is 



g N (p',z) = 



1 



JV-l 



N 



1 ^ 

n=0 



i b {n 7 z)AI(p',n,z), (23) 



where, anticipating the need to approximate the DC-block 
used in continuous- wave pseudothermal ghost imaging to 
obtain a high-contrast image Q , we have chosen to cross- 
correlate 



N-l 



i b (n, z) = i b (n, z) - ^- ^ i b { 



(24) 



with 



AI(p', n, z) = I{p', n, z) - (I(p\ n, z)), (25) 

for which the ensemble average needed in the second term 
can be computed because we know the statistics of the 
spatial-pattern sequence n)} being applied to the 

SLM. 



B. Spatial Resolution 

The computational ghost imager's spatial resolution is 
found from the ensemble average of Eq. ([23)) . We will fo- 
cus our attention on the ranges at which there are target 
components, i.e., z€{z k :l<k<K}. Furthermore, 
we shall assume that the ghost imager's pulse duration is 
short enough to resolve all these ranges; for our p(t) this 
is guaranteed if cT p /2 < mmi<k<K-i (zfc+i — z k ). Thus 
we need to evaluate (Gn(p' , z k )) , which simplifies to 

(Gn(p',zic)) = (i b (n,z k )I(p',n,z k )) 

- (i b (n,z k ))(I(p',n,z k )), (26) 

for 1 < k < K, where we have employed the {£(p, n)} 
being a sequence of statistically independent, identically 
distributed, spatial patterns. 

The mean values that appear in the preceding expres- 
sion are easily computed. For the computed reference we 
have 



(I(p',n,z k )) 



dpi J dp 2 (f(pi,n)£{p 2 ,n)) 



k 2 e 



ik [\ Pl \ 2 -\p 2 \ 2 -2p'-( Pl -p 2 )]/2z k 



(2nz k y 



-2|pT/< 



(27) 



where a Zk = \oz k /irpo is the intensity radius of 
\E Zk (p', t)\ 2 , and the approximation follows from our as- 
sumptions that po <C &o and koaopo/2z k <C 1. For the 
bucket detector we find 



(i b (n,z k )) = qi]N T 



J dr g(nT s + 2z k 



r) 



dp b Al(p b )(\E d (p b , 



2g V N T A b 



dp'T Zk (p')e 



-2|p'l 2 /^ 



(28) 



,(29) 



where po ^ min(a , psh 
tion, A b = J dp b A 2 b {p b ) 



/Ws) justifies the approxima- 
irr 2 1 2 is the area of the bucket 
detector's photosensitive region, and T' = l/ Jdt\p(t)\ 4 
is the effective duration of the matched filter's response 
to |p(£)| 2 . The exponential term appearing in Eq. (f2"9")l 
comes from the average intensity pattern in the z k plane, 
as seen in Eq. (f27|) . Going forward, we shall assume that 
the target lies well within the center of this average in- 
tensity pattern, which reduces Eq. (|2"T| to 



(I(p',n,z k )) 



on the range-Zfc target, and Eq. (|29|) to 



(i b (n, z k )) 



2qi]N T A b 



T'z 2 



dpT Zk {p). 



(30) 



(31) 



To complete our evaluation of the average range-z^ 
ghost image, (Gn(p' ', z k )) , we need derive an expression 
for (i b (n, z k )I(p' ,n, z k )) . Its derivation is more involved 
than what sufficed above for (I(p' ,n, z)) and (i b (n, z k )) , 
but, because it parallels similar analysis from [13j . we 
shall merely describe the procedure and present the re- 
sult. First, we backpropagate the bucket-detector mea- 
surements to the z — source plane. Next, we exploit the 
statistical independence of the sequence of spatial pat- 
terns, the turbulence present on the source-to-target and 
target-to-receiver paths, and the target's field-reflection 
coefficient. We then employ our far-field propagation 
conditions, the correlation function of the target's field- 
reflection coefficient, and the mutual coherence functions 
for the turbulence fluctuations incurred on the two paths. 
Finally, to evaluate a fourth-order moment of the spatial 
patterns £(p, n) we employ Gaussian moment-factoring, 
and thus express the fourth-order moment in terms of 
second-order moments given in Eqs. © and ([3]). The 
average ghost image term in Eq. ([26]) is then found to be 



{G N (p',z k )) = 



AqrjNTA b p 



dpT Zk {p)- 



-\p' -P\ 2 /apl. 



(32) 

with p Zk = XoZk/wao being the range-z^ coherence length 
of the transmitted spatial patterns {£(p, rt)}, and a = 
1 + agWs/2p| the resolution-degradation factor imposed 
by the turbulence present in the source-to-range-Zfc path. 
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Equation (|32|) shows that the average computational 
ghost image for the target component at range z k is 
T Zk (p) convolved with a Gaussian point-spread function 
(PSF) of width y/ap Zk . In the absence of turbulence 
a = 1, so that we get the spatial resolution previ- 
ously found for single-range, continuous-wave operation 
in [8(, viz., the speckle coherence length. In the pres- 
ence of turbulence in the source-to-targct path we have 
a > 1, and our result coincides with that for single-range, 
continuous-wave operation through turbulence given in 
[l3j . The key points to be gleaned from Eq. (j52")l are: (1) 
turbulence in the source-to-target path does not signifi- 
cantly degrade spatial resolution until its spherical-wave 
coherence length in the source plane, ps/VWs, becomes 
smaller than the source's intensity radius, ao; (2) the 
spatial incoherence of the rough-surfaced target leads to 
turbulence in the target-to-receiver path having no ef- 
fect on the average ghost image; (3) computational ghost 
imaging does not suffer from the turbulence that a pseu- 
dothermal configuration would were its propagation path 
through turbulent air; and (4) our correlating %{n,z) 
with AI(p',n,z) has eliminated the background term 
from (Gn(p' , z)) — which would have limited image con- 
trast had we instead correlated ib{n, z) with I{p' , n, z) — 
resulting in a high-contrast image. 

A final point to be made here concerns depth of fo- 
cus and coherence propagation. Our far-field assump- 
tion for the ghost imager, k^a^p^jlzk "C 1, is an ana- 
lytic tool that simplifies the propagation of the ensemble- 
averaged correlation function in Eq. ([3]) from the source 
to range Zk- However, it does not imply that specific 
field patterns can be calculated with a far-field approxi- 
mation, i.e., a spatial Fourier transform, so that the in- 
tensity patterns at different ranges differ only by coor- 
dinate and amplitude scaling. Because po <C ao, there 
is a significant region wherein koaopo/2zk <C 1, but we 
are in the near field of the source's intensity diameter, 
viz., fcoa§/2zfc 3> 1. In this regime, it is important that 
the reference used to form the range- Zk ghost image be 
computed specifically for that range. This is because 
computational ghost-image spatial resolution, in the ab- 
sence of turbulence, is known to degrade by a factor of 
y/l + (Sz I Zk) 2 (koaQ/ '4zfe) 2 , where Sz = z — Zk, when the 
reference used was computed for range z instead of range 
Zk [H- So, deep in the near field a small range mismatch, 
\Sz\/zk <C 1, between the reference and the target can 
substantially degrade the spatial resolution. The same 
degradation applies to a pseudothcrmal ghost imager, in 
which a physical reference arm is used, hence in that case 
a different reference measurement must be made for each 
target range of interest if the range spread exceeds that 
imager's depth of focus. 



C. Signal-to-Noise Ratio 



its variance, 



SNR eN (p',z fe ) 



(G 2 N (p',z k ))-(g N (p',Zk)y 



(33) 



To simplify the derivation we make several additional 
assumptions. First, that the ghost imager resolves 
all significant detail in the target, allowing us to use 

fdpT Zk (p)e~ lp '~ pl2/af>l * ~ Trap 2 Zk T Zk (p'). With this as- 
sumption, the square of Eq. (|32p becomes 



(G N {p',z k )f 



Aqi]N T A b pl 



(34) 



This leaves the tedious task of finding (G%(p', z k )), which 
involves: eighth-order and sixth-order moments of the 
{£(/>, n)}; propagation through turbulence to the target 
element at range Zk, including a fourth- moment evalua- 
tion for the incurred turbulence fluctuations; target re- 
flection, including a fourth-moment evaluation for the 
field-reflection coefficient; propagation back to the re- 
ceiver, requiring another turbulence-fluctuation fourth 
moment evaluation; and photodetection, with is accom- 
panying shot noise and matched filtering. The statistical 
independence of the sequence of spatial patterns, the tur- 
bulence fluctuations on each path, and the target's field- 
reflection coefficient is a considerable help in completing 
this evaluation. So too are the Gaussian distributions of 
the {£(p, n)} and T Zk (p / ), which permit high-order mo- 
ments to be found from the second moments we have 
presented earlier via Gaussian moment-factoring. 

The most vexing difficulty in the calculation turns 
out to come from the fourth moments of the turbu- 
lence fluctuations. Here we shall follow the lead provided 
in [f3j], where the corresponding SNR for single-range, 
continuous-wave operation was derived. Our work, how- 
ever, will differ from that in [I3| in that it will account for 
pulsed operation in which a single pulse's duration, T p is 
much shorter than the atmospheric coherence time, but 
NT S , the time duration of the iV-pulse sequence used to 
form the ghost image, will be taken to span 1 < N c < N 
atmospheric coherence times. 

Following [l3| we take the logamplitude fluctuations to 
be Gaussian distributed, and assume that the turbulence 
is of weak-to-moderate strength — or sufficiently concen- 
trated near the target — that we can both ignore the coor- 
dinate dependence of the turbulence at the transmitting 
and receiving planes, and assume that the logamplitude 
coherence length at the range-Zfc target is larger than p Zk . 
Under these conditions, the relevant turbulence-fourth 
moments for the source-to-target and target-to-detector 

paths reduce tO e ^%K S [(n-n')T s \ &nd e 4a 2 D K D [(n-n')T s ] 

spectively, for the correlation between turbulence affect- 
ing pulses n and n'. Here, for m — S, D, K m (r) is the 
logamplitude fluctuation's normalized (K m (0) = 1) co- 
variance function, and 



We will now evaluate the ghost imager's signal-to-noisc 
ratio, defined to be ratio of Gn(p' , z&)' s squared mean to 



0.562 k 



7/6 



dzCl {z) 



z{z k - z) 



Zk 



5/6 



(35) 
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is its Rytov-approximation variance [24[ . To account for 
temporal averaging of the turbulence, we define 



N-l 



2 s K s (nT s )+4<T 2 D K D (nT s ) _ 1 



n=0 



(36) 



For N 2> 1 in the assumed weak-to-moderate 
turbulence — for which er^ < 0.1 — we will have 7 = 
(e 4cr s +4cr o - \)/N <C 1 when the turbulence decorrc- 
lates pulse-to-pulse on both paths, so that Ks{nT s ) = 
Ku(nT s ) = (5„o- On the other hand, when the turbu- 
lence is frozen across all N pulses, so that Kd(tiT s ) = 
Ks(nT s ) = 1 for all 71 of interest, we get 7's other asymp- 
tote, 7 = e 4o " s+4CTD — 1. Finally, to simplify the analysis, 
we will assume the transmitter pulse is flat-topped, viz., 
p(t) = for -T p /2 < t < T p /2 0. 

The preceding tools can now be employed to show that 



SNR e „(p',z fc ) = 



A Z S + A Z TZ + A Z V + A J 



(37) 



with the terms that appear in the noise denominator be- 
ing 



A 2 S 



A TZ 



A 2 V 



AT 



2 , „ 1 + 2 7 (1 + 0) 



1 + 2/3 

Zfc^ZS+fc 



2A b irr]N T Np 4 Zk ' 

A* h (f*(«B+*B> - 1-7) 



(38) 
(39) 
(40) 

(41) 



where A Zk = J dp T Zk (p) and A' Zk = J dp 7? (p) are two 
measures of the target's area, and j3 = r 2 /a^ measures 
the bucket detector's area relative to the source's inten- 
sity area. 

These noise terms account for the following phenom- 
ena. The A 2 S term is the noise contribution from the 
source-produced on-target speckle patterns. That it de- 
creases inversely with the number of pulses, N, used to 
form the ghost image is indicative of the need to use 
many different illumination patterns to form a ghost im- 
age. This noise term grows linearly with the ocA' /ftp 2 
number of spatial-resolution cells on the target, and is ex- 
acerbated by the presence of turbulence, i.e., c|+er|, > 0. 

The A 2 R term in the noise denominator arises from the 
random fading that results from the rough-surface target 
reflection combined the effects of atmospheric turbulence. 
When NT S spans N c ^> 1 turbulence coherence times, we 
get 7 — > 0, indicating that the turbulence contribution 
to this noise term vanishes. However, because we have 
not allowed our rough-surfaced target to decorrclate over 
this measurement interval, its contribution to this noise 
term can only be reduced by aperture averaging, i.e., by 



increasing the bucket detector's area so as to capture and 
average an increasing collection of uncorrelated target 
speckles. 

The A 2 T> noise term is the bucket detector's shot-noise 
contribution, which is inversely proportional to the aver- 
age number of detected target-return photons. For fixed 
optics and a given target, this term is only decreased by 
increasing the average transmitted photon-number per 
pulse, Nt, or the number of pulses, N, used to form the 
ghost image. 

The preceding noise terms were present — albeit with 
somewhat different scaling factors — in the SNR expres- 
sion for continuous-time ghost imaging [l3j , but the A 2 J 7 
noise-denominator term in Eq. (|37|) is a heretofore unen- 
countcred consequence of time-varying turbulence dur- 
ing the measurement interval. When the turbulence is 
frozen across all N pulses, this term disappears, but when 
the turbulence changes during the measurement it cre- 
ates randomness in the on-target illumination patterns 
that changes from pulse to pulse. This new randomness 
is suppressed by increasing the number of pulses used 
to form the ghost image. To appreciate the impact of 
this new term, consider the worst-case scenario, in which 
the turbulence decorrclates pulse-to-pulse, the new term 
dominates the noise denominator, and N is to be cho- 
sen to achieve a desired SNR while still under the sway 
of this term. Here we find that the necessary N value 
is proportional to the square of the number of spatial- 
resolution cells on the target. In contrast, were there no 
time variation to the turbulence — so that the fourth noise 
term vanished — then N would only need to be propor- 
tional to the number of on-target spatial-resolution cells 
to achieve the desired SNR when the first and/or third 
noise terms dominate the second. Because increasing N 
with T s fixed increases the image acquisition time, the 
fourth noise term can have a significant adverse impact 
on that acquisition time. 

Having described the various noise contributions to the 
ghost-image SNR, let us conclude this section by high- 
lighting two physically-important asymptotic forms of 
Eq. (f3"T)l . When the number of pulses increases with- 
out bound, the ghost-image SNR reaches a finite limit, 
called the saturation SNR, which is given by 



SNR Sjv , sat (p',z fc ) = 1 + 2/3. 



(42) 



Note that it is due solely to the time-independent tar- 
get speckle [26|, and is independent of the target's range 
and reflectivity, although the Nt and N values needed to 
reach the saturation regime do depend on 7~ Zk {p')- This 
saturation SNR can only be increased by increasing /3, 
which means either increasing the size of the bucket de- 
tector, or decreasing the size of the source, with the lat- 
ter entailing a degradation of the ghost imager's spatial- 
resolution capability. 

When the average number of detected target-return 
photons is so low that shot noise dominates all other 
fluctuations in the ghost image, we get the shot-noise- 
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limited SNR 

SNR gjv , shot (p ,Zk) = j 2 2 • ( 43 ) 

Unlike the case for its saturation signal-to-noise ratio, 
the ghost imager's shot-noise limited SNR does depend 
on the target's range and reflectivity. 

IV. LASER RADAR 
A. Image Formation 



the discussion at the end of Sec IIIIBI about this issue 
for the ghost imager. When k r]/2z <C 1 for all target 
ranges of interest, the laser radar can be focused at in- 
finity with no loss of spatial resolution at any of those 
ranges. However, when korf/2z ^> 1 for target ranges of 
interest, then the radar receiver must operate within the 
depth of focus for the range to be imaged in order to pre- 
vent loss of spatial resolution. The rest of our treatment 
of the laser radar will assume that the far-field condition, 
k r^/2z -C 1, holds, so that all target ranges of interest 
will be in focus, but this need not always be the case in 
operational scenarios of interest. 



i p (p' p7 n,z)= drg(nT s +2z/c-T) 



The -/Vp-pixel laser radar image for range z is formed 
from the pixel-wise outputs { i p {p' p , n, z) : 1 < p < 
N p , < n < N - 1} of a CCD— on which the target- 
reflected light has been focused by a lens — that have been 
matched-filtered for range z. These photocurrents take 
the form 



(44) 



where: A p (p' i ) is the real- valued field-transmission pupil 
function defining the photosensitive region of the pth 
image-plane pixel, which we take to be centered at p\ = 
—p'p to compensate for image inversion; E p {p' il r) is the 
complex envelope of the light impinging on the CCD 
plane at point p[ and time r; and Ai p (r) is the pho- 
tocurrent shot noise for pixel p. 

In practice the lens' focal length will cast a minified 
image on the CCD, but, for convenience, we will assume 
it is chosen to realize 1:1 imaging for range Z£, which is 
taken to be the center of the range-interval of interest. It 
follows that 



E p (p',t) = dpAt{p)E d {p,t) 



x e 



-zk \p\ 2 /z e foe 

i2nz£ 



, (45) 



where Ae(p) is the lens' real- valued field-transmission 
pupil function, Ed{p,t) is the pupil-plane target-return 
field given by Eq. (|19|) . and we have neglected an unim- 
portant absolute phase factor as well as the propagation 
delay within the radar receiver. To facilitate analytic 
comparison with the ghost imager, we will take the lens' 
pupil function to be Ai{p) = exp(— \p\ 2 /rj). The laser 
radar then produces its range- z image by pixel-wise av- 
eraging the photocurrents from the N pulses, 



N-l 



£n(p' p ,z) = — i P (Pp,n,z) 



(46) 



Before turning to our spatial-resolution analysis, there 
is an important point to make about depth of focus, cf. 



B. Spatial Resolution 

The laser radar's spatial resolution is found from its 
average image, just as was done for the ghost imager. As 
we did in Sec. lIIIBl we will only consider ranges {z k : 1 < 
k < K } that contain target components, and take the 
pulse duration, T p , to be short enough to resolve them. 
We then get 

(£ N (p'z k )) = q-q dr g(nT s +2z k /c- t) 



x jdp'AKtimM^r). (47) 

Backpropagating E p (p' , r) to the source using Eqs. P5|) . 
([T^| . and JS]), evaluating the resulting source and target 
second-moments that appear, and then performing all 
integrations but those over the target and image planes, 
we obtain 



(C N {p' p ,z k )) 



2q?]N T A e 



x dp'iAlip'i)- 



dp'T Zk (p')e-^'^< 

p'+p'i\ 2 / a 'pJi 



(48) 



In this expression: w Zh — XoZk/trwo is the radar beam's 
intensity radius at range z k , which we will assume is suf- 

ficiently large that % k (p')e~ 2|p ' |2/ <= » T Zk (pJ, A t = 
JdpAj(p) = Ttrf/2 is the lens' effective area; p' = 

Xo z k/ '\f2itri is the diffraction- limited (no-turbulence) 
spatial resolution for an image-plane point detector; and 
a' = 1 + t^Wd/Pd 1S a resolution-degradation factor 
caused by turbulence. Turbulence begins to impair the 
laser radar's spatial resolution when its receiver-plane co- 
herence length, pi) j \JWd , becomes comparable to the 
receiver's lens radius, rg. Hence, when the turbulence 
on the return path is either sufficiently weak, or concen- 
trated near the target, we get a 1 as 1 and there is no 
loss of resolution. However, for sufficiently strong turbu- 
lence, or turbulence concentrated near the receiver's lens, 
significant resolution degradation is incurred. 

The computational ghost imager can calculate its ref- 
erence field for every p' of interest, but the laser radar 
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must use pixels of finite size in order to collect any optical 
power. To obtain a closed-form result from the pixel in- 
tegration in Eq. (|4"5)l we will use Ap(p') = e~' p +p p l T 'v 
and arrive at our final expression for the laser radar's 
average image, 



(C N (p' p ,z)) 



2q7]N T AiA F 



dp'% k (p')- 



3 -|p'-Pp| 2 /(/3'+a')p4. 



tt(/3' + a')p[ 



(49) 



where A p = Jdp' A p (p r ) = Trr p /2 is the pixel's effective 

area, and f3' = r p /2p z 2 k . Equation (|4"9")l quantifies the 
laser radar's loss of spatial resolution — its degradation 
from p' Zk , the diffraction-limited, point-detector value — 
when the lens area exceeds a turbulence coherence area 
and/or the pixel area exceeds the diffraction- limited spot 
size. 



C. Signal-to-Noise Ratio 

The laser radar's SNR for the pth pixel's range-z^ value 
is the ratio of £n{p p , Zfc)'s squared mean to its variance, 



SNR £jv ( P ;,z fe ) = 



(£ 2 N (p pl z)) - (C N (p p ,z)Y 



(50) 



Our evaluation of this SNR mirrors what we did for the 
ghost imager. So, as we did there, we shall assume that 
all significant target detail is resolved by the imager, and 
that the turbulence is weak enough that we can ignore its 
fluctuation-terms' coordinate dependence in the planes 
of the transmitter, target, and receiver. With the first 
assumption, Eq. (|49|) yields 



(C N {p' pl z k )) 2 = 



2qT 1 N T AiA I 



TTZ 2 k W 2 z T p 



Tz k (p' p ) 



(51) 



As was the case for the ghost imager, the primary diffi- 
culty encountered in SNR evaluation is finding the im- 
age's second moment. The laser radar's second-moment 
calculation, however, is substantially simpler than that 
for the ghost imager in that the laser radar case only 
requires turbulence and target fourth-moments and the 
shot noise's second moment. Consequently, the evalua- 
tion proceeds as follows. 

First, we use the statistical independence of the tur- 
bulence, target, and shot noise to separate averages in- 
volving these three randomness contributors. Then we 
employ Gaussian moment-factoring to reduce the tar- 
get's fourth moment into a sum of products of second 
moments, and our assumption about the coordinate- 
independence of the turbulence fluctuations to evalu- 
ate that fourth moment as we did for the ghost imager. 
Next, we perform the resulting multi-dimensional Fourier 



transforms of the Gaussian functions that arise from far- 
field optical propagation in conjunction with the Gaus- 
sian pupil functions we have assumed. Finally, under our 
assumption that the laser radar resolves all significant 
target detail, we are left with 



SNR £N ( P ;,z fe ) 



r Zk (p' P )( 



l + (2+p"b > 
1+13' > 



2r)N T NA t A p 



(52) 



where 7 is the time-averaged turbulence factor from 
Eq. (HO). 

The noise-denominator terms in the preceding SNR 
formula have the following physical interpretations. The 
first term is due to the time-independent target speckle, 
exacerbated, to some degree, by the turbulence-induced 
scintillation on the transmitter-to-target and target-to- 
receiver paths. The second term is due to shot noise. 
Thus when NtN — the average number of photons trans- 
mitted over all TV pulses — is sufficiently high, the SNR 
reaches a finite maximum value, namely the saturation 
signal-to-noise ratio, given by 



SNR£ JV , sat (p' , Zk) = 1 + /3', 



(53) 



which depends only on r 2 /2p 2 k , the number of 
diffraction-limited spots within the pixel area (26[. In- 
creasing (3' , e.g., by increasing the pixel size, will increase 
the saturation value, but doing so degrades the radar's 
spatial resolution; sec Eq. (|4"9")l . Conversely, when the 
SNR is much lower than its saturation limit, it takes the 
shot-noise limited form 



SNR 



C n ,shot 



(Pp,2k) 



2r,N T NA,A p T Zk (p' p ) 



(54) 



V. PERFORMANCE COMPARISON 

Having completed spatial resolution and SNR analyses 
for both the computational ghost imager and the laser 
radar, we are ready to compare their capabilities. Before 
proceeding, two points deserve note. The first concerns 
ensuring that our spatial resolution and SNR compar- 
isons are fair. Toward that end we will take wq = po, so 
that the laser radar's on-target intensity pattern matches 
the ghost imager's average on-target intensity pattern. 
We will also assume that ri = a , because: (I) n and 
ao correspond to intensity radii (r£ for the the intensity 
transmission of the radar receiver, and ao for the intensity 
transmission of the ghost imager's exit optics); and (2) 
ri and ao determine the diffraction-limited spatial reso- 
lutions of these two imagers. In addition, we we will take 
Nt gi Nqi = NtluNlr, i-e., the product of the average 
number of transmitted photons per pulse and the num- 
ber of pulses employed for the ghost imager (GI) and the 
laser radar (LR) must be equal, but we will not require 
Nt gi = Nt lr and -/Vgi = -/Vlr, as originally indicated 
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in Sec. ITT] Thus we are constraining both systems to use 
the same average number of transmitted photons to form 
their images, but we are affording the laser radar the op- 
portunity to exploit its ability to form an image from a 
single pulse, in contrast to the ghost imager's fundamen- 
tal requirement for averaging the returns from multiple 
pulses. 

Our second pre-comparison note is cautionary. Our 
analysis has employed Gaussian functions for partially- 
coherent and coherent optical beams, and for photode- 
tector pupil- functions, etc., in order to obtain closed- 
form expressions for our systems' spatial resolutions and 
SNRs. Although the use of Gaussian functions should 
yield the correct dependence of these performance met- 
rics on system parameters, the constant factors that ap- 
pear in our results would be different had we numerically 
evaluated the spatial resolutions and SNRs for circular- 
pupil optics. 



A. Spatial-Resolution Comparison 

We have shown that the ghost imager and laser radar's 
average images for the target component at range Zk con- 
volve T Zk (p') with Gaussian point-spread functions of the 
form P(p') = e~' p ' / r ™/7rrf cs , where 



'Otp z 



+ a' p 



2fc' 



ghost imager 
laser radar, 



(55) 



with p Zk and p' Zk being the diffraction-limited resolutions, 
a > I and a' > 1 accounting for turbulence-induced 
resolution loss, and /?' > for resolution lost because of 
finite pixel size. 

The diffraction- limited resolutions are p Zk = XoZk/nao 
and p' Zk = AoZfe/v^ 7rao, where we have used re = ao 
in the latter, showing that the laser radar's diffraction- 
limited performance is -s/2-times better than that of the 
computational ghost imager. The numerical value of this 
resolution advantage depends on our use of Gaussian 
functions, but the fact that the laser radar's diffraction- 
limited spatial resolution is better than that of the com- 
putational ghost imager does not, so long as the laser 
radar's receiver pupil is the same as the ghost imager's 
transmitter pupil. That this is so follows from the ghost 
imager's PSF arising from convolution of the on-target 
average intensity pattern of its transmitter with the cor- 
responding pattern from the computed reference, while 
the laser radar's PSF (in the 1:1 imaging setup we have 
assumed) corresponds to just one of those patterns. 

As noted in Scc. lIVl the laser radar must employ finite- 
sized pixels, and Eq. (|52p shows that — other system pa- 
rameters being held constant — forcing /?' <C 1 will push 
the system into its shot-noise limited regime as the pixel 
decreased. Thus a prudent compromise might 
be to size the pixels to satisfy /?' = 1, so that the laser 
radar's spatial resolution, in the absence of turbulence, 
only suffers a v^-factor increase. This choice of pixel 



size, however, wipes out the laser radar's spatial resolu- 
tion advantage over the computational ghost imager, i.e., 
the no-turbulence resolution of the former then exactly 
matches the diffraction-limited resolution of the latter. 

When turbulence is strong enough to dominate the 
spatial resolution of both imagers, we get 



''re 



XoZkVWs/nVZps, ghost imager 
XoZk^/WB /fty^ Pd, laser radar, 



(56) 



where we have continued our use of re = clq. For a sit- 
uation in which ps/VWs = Pd/VWd — such as when 
the laser radar's receiver is co-located with the ghost- 
imagcr's transmitter — we find the turbulence-limited res- 
olution to be equal for both systems. More gener- 
ally, either system could have more favorable spatial- 
resolution behavior in the presence of turbulence because 
different turbulence-strength profiles could exist on the 
transmitter-to-target and target-to-receiver paths, and 
the ghost imager's spatial resolution is only sensitive to 
turbulence on the transmitter-to-target path while the 
laser radar's spatial resolution is only sensitive to the 
turbulence on the target-to-receiver path. 



B. Signal-to-Noise Ratio Comparison 

For our SNR comparison we shall consider the two 
imagers' saturation SNRs and their shot-noise limited 
SNRs. Throughout we will assume that any turbulence 
that might be present has no effect on spatial resolution, 
so that a = a' = 1, and that both imagers have sufficient 
resolution to resolve all significant detail in T Zk (p')- 

Equations ([42]) and (|53|) specify the saturation SNRs 
for the ghost imager and the laser radar, and give us 



SNR £iv , sat (p;,z fc ) _ \ + fj' 
SNR Sjv , sat (p^,z fc ) ~ 1 + 2/3' 



(57) 



Setting rf — ao, to make the laser radar's receiver pupil 
be the same size as the ghost imager's transmitter pupil, 
and /3' = 1, so that these systems have the same spa- 
tial resolution, the preceding ratio of their saturation 
SNRs becomes 2/(1 + j3). Consequently, the compu- 
tational ghost-imager's saturation SNR will greatly ex- 
ceed that of the laser radar when j3 = r\ja\ 3> 1, i.e., 
when the bucket detector's area is much larger than the 
source area. There is a simple physical explanation for 
this behavior: the saturation SNR is due to the time- 
independent speckle created by reflection from the rough- 
surfaced target. These speckles arc ^ao in radius in the 
receiver's pupil plane, so when (3 1, the bucket de- 
tector is averaging over many statistically independent 
speckles, driving up the ghost imager's saturation SNR. 

There is, however, an intrinsic unfairness in the pre- 
ceding favorable view of the ghost imager's saturation 
SNR, because by fixing re = ao but letting j3 1 , we are 
allowing the bucket detector to have a much larger re- 
ceiving aperture than the laser radar. Suppose, instead, 
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that we constrain the systems to have the same size re- 
ceiving aperture, so that = r>, but maintain /3 > 1. 
In that case the laser radar's diffraction-limited spatial 
resolution would be better than that of the ghost imager 
by a factor of r^/ao 3> 1 but its SNR would be worse. 
However, the laser radar could then use /3' = 2r\ja^ ^> 1 
and: (1) have its no-turbulence spatial resolution match 
the diffraction limit of the ghost imager; and (2) have 
its saturation SNR equal that of the ghost imager. In 
this case, the time-independent speckle is being averaged 
in the laser radar's image plane, because the large pixel 
comprises a great many independent speckle lobes. 

Finally, if we make all optics the same size, ao = r& = 
m, and choose the pixel size for equal spatial resolu- 
tions, then f) = f3' = 1 and we find SNRc N , S at(p' p , zt) = 
2SNR< ?N , sat ( / ^,z fc )/3. 

To summarize what we have seen so far concerning 
saturation-SNR behavior, fair comparisons between the 
ghost imager and the laser radar indicate that neither 
system will enjoy a significant advantage in this regard. 
There is, however, a fine point to be considered concern- 
ing atmospheric turbulence. We have assumed the rele- 
vant turbulence coherence lengths to be long enough that 
our imagers' spatial resolutions do not lose any resolution 
arising from propagation over turbulent paths. Neverthe- 
less, the target-speckle terms in the noise denominators 
of Eqs. (|57|) and (f52"|) both contain the time-averaged scin- 
tillation factor 7 When a sufficient number of pulses are 
averaged — as has been assumed in the saturation SNR 
formulas — we get 7 — > 0. For the ghost imager, averaging 
the returns from a large number of pulses is an intrinsic 
requirement for image formation, but a laser radar can 
form its image with a single pulse, which case we get 

SNR £l ( P ;,z fc )^ TT i±^-, (58) 

for N Tlr sufficiently high. For /3' = 2rf/a§ > 1, this 
SNR can be substantially worse than the ghost imager's 
saturation SNR, because 7 > e Aa " s in this single-pulse 
case. The laser radar will then need to employ more 
than a single pulse to approach its full saturation SNR — 
and to match that of the computational ghost imager at 
the same spatial resolution — because of scintillation. 

Now let us turn to SNR behavior when neither system's 
Nt and N values are sufficient to reach SNR saturation 
by comparing their shot-noise limited performance. Here, 
using ri = ao, Eqs. P5|) and ([53)1 give us 

SNR £jv , shot (p;,z fc ) = A Zk /3' 
SNR Whot (p^ fc ) 2K(% h P% h (p> p y 1 ' 

For a fair comparison we again set /3 = /?' = 1, to give the 
systems the same pupil sizes and spatial resolutions. The 
preceding SNR ratio then reduces to A Zk /2TTp 2 Zf T Zk {p' p )- 
For a range-Zfc target whose intensity-reflection coeffi- 
cient has limited spatial variation, this number is ap- 
proximately its number of spatial-resolution cells, mak- 
ing the laser radar's shot-noise limited SNR far superior 



to that of the computational ghost imager. Increasing 
the size of the ghost imager's bucket detector can over- 
come this SNR disadvantage by making ^ > 1. As in 
our saturation-SNR comparison, however, we should then 
allow the laser radar to increase its lens size to equal 
the ghost imager's new value, while operating with 
f3 = r^/af) 3> 1, to regain its shot- noise limited SNR 
advantage. Unlike what we found for the laser radar's 
saturation SNR, its shot-noise limited SNR is the same 
at all -/Vt lr -/Vlr values, i.e., Eq. (|59|) applies even when 
the laser radar forms its image from a single pulse. 



VI. CONCLUSIONS 

Computational ghost imaging is, in many respects, a 
dual of floodlight-illumination laser radar. The com- 
putational ghost imager's system complexity is in its 
transmitter, whose sequence of SLM patterns creates the 
structured illumination that provides the imager's spatial 
resolution. The laser radar's complexity lies in its CCD 
receiver, which provides its spatial resolution. Conse- 
quently, the size of the ghost imager's transmitter pupil 
sets its diffraction-limited spatial resolution, whereas the 
laser radar's no-turbulence spatial resolution is set by the 
size of its receiver pupil in conjunction with that of its 
CCD pixels. Thus only turbulence on the transmitter- 
to-target path can impair the ghost imager's spatial res- 
olution, while the laser radar's spatial resolution is only 
impacted by turbulence on the target-to-receiver path. 
Therefore, ghost imaging and laser radar systems that 
are designed to have equal spatial resolutions in the ab- 
sence of turbulence could have significantly different per- 
formance in a bistatic configuration, in which their trans- 
mitters and and receivers are not co-located, so that sig- 
nificantly different turbulence distributions are encoun- 
tered on these two paths. In such situations cither one 
could offer the better spatial-resolution performance, de- 
pending on which path had its turbulence concentrated 
near the resolution-controlling pupil. Aside from this tur- 
bulence issue, our analysis indicates that a fair compar- 
ison between the computational ghost imager and the 
floodlight-illumination laser radar shows them to have 
equal spatial resolutions, except for the following two 
caveats: (1) the laser radar can form its image from a 
single pulse, making it far better for imaging moving tar- 
gets; and (2) the computational ghost imager has infinite 
depth of focus, whereas the laser radar will not for ranges 
satisfying fcor?/2z ^> 1. 

Both the ghost imager and the laser radar have signal- 
to- noise ratios that are shot- noise limited at low-iVr, low- 
TV values and saturate at high- iVr, high- TV values. When 
their optics are sized for equal spatial resolutions, with 
Nt Q i = Nt lr and Nqi = TVlr, there is little difference 
in their saturation SNRs. This contrasts strongly with 
their shot-noise limited SNR behavior, under these con- 
ditions, in which the laser radar outperforms the ghost 
imager by a factor approximately equal to the number 
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of spatial-resolution cells on the target. As a result, we 
can expect that the ghost imager will require significantly 
more time to achieve a desired SNR when operating in 
this regime. This key disadvantage for correlation-based 
ghost imaging could be mitigated to some degree, how- 
ever, by the use of compressed-sensing techniques, which 
enable many fewer pulses to suffice for ghost-image for- 
mation. Recent work has demonstrated this possibility in 
table-top ghost imaging done in reflection |27| , indicating 
its likely feasibility for standoff-sensing applications. 

What then are the possible advantages of ghost imag- 
ing in comparison with laser radar? The principal 
such advantage identified by our analysis accrues in 
bistatic configurations wherein, for operational reasons, 
the transmitter must be located in a region of weak 
turbulence but the receiver necessarily is in a strongly- 
turbulent region. Beyond that, however, there are some 
technological possibilities. The ghost imager only re- 
quires a single-pixel detector, whereas the laser radar 
needs a detector array. For wavelength regions in which 
high-performance single-pixel detectors arc available but 



similar-quality detector arrays arc not, ghost imagers 
would provide active- imaging capability that laser radars 
could not. A related technological advantage arises for 
ghost imaging in multi-static configurations, in which a 
network of simple, small, single-pixel detectors view a 
target region that is floodlit by a single, large-aperture, 
structured-illumination transmitter. Individual images 
could be formed from each detector's outputs to capture 
multiple views of the target, and allow for more averag- 
ing of the target-induced speckle. A corresponding multi- 
static laser radar would require high-resolution CCDs at 
each receiver location, making it more complicated and 
more expensive than the ghost imager. 



Acknowledgments 

This work was supported by the DARPA Information 
in a Photon Program under U.S. Army Research Office 
Grant No. W911NF-10-1-0404. 



[1] B. I. Erkmen and J. H. Shapiro, Adv. Opt. Photon. 2, 
405 (2010). 

[2] T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. 

Sergienko, Phys. Rev. A 52, 3429 (1995). 
[3] G. Scarcelli, V. Berardi, and Y. Shih, Phys. Rev. Lett. 

96, 063602 (2006). 
[4] F. Ferri, D. Magatti, A. Gatti, M. Bache, E. Brambilla, 

and L. A. Lugiato, Phys. Rev. Lett. 94, 183602 (2005). 
[5] B. I. Erkmen and J. H. Shapiro, Phys. Rev. A 77, 043809 

(2008) . 

[6] B. I. Erkmen and J. H. Shapiro, Phys. Rev. A 79, 023833 

(2009) . 

[7] J. H. Shapiro and R. W. Boyd, Quantum Inf. Process. 

11, 949 (2012). 
[8] J. H. Shapiro, Phys. Rev. A 78, 061802 (2008). 
[9] D. Venkatraman, N. D. Hardy, F. N. C. Wong, and J. H. 
Shapiro, Opt. Lett. 36, 3684 (2011). 
[10] O. Katz, Y. Bromberg, and Y. Silberberg, Appl. Phys. 

Lett. 95, 131110 (2009). 
[11] R. E. Meyers, K. S. Deacon, and Y. Shih, Phys. Rev. A 

77, 041801 (2008). 
[12] C. Zhao, W. Gong, M. Chen, E. Li, H. Wang, W. Xu, 

and S. Han, Appl. Phys. Lett. 101, 141123 (2012). 
[13] N. D. Hardy and J. H. Shapiro, Phys. Rev. A 84, 063824 
(2011). 

[14] N. D. Hardy and J. H. Shapiro, Proc. SPIE 7815, 78150L 

(2010) . 

[15] B. I. Erkmen, J. Opt. Soc. Am. A 29, 782 (2012). 

[16] J. Cheng, Opt. Express 17, 7916 (2009). 

[17] Ghost imagers intrinsically require multiple pulses to 
form an image, whereas a laser radar can provide an 
image from a single pulse. For now we will allow them 
to use the same number of pulses, but we will return to 
this distinction, in Sec.[V] after we have derived the two 
systems' spatial resolution and SNR behaviors. 



[18] M. I. Skolnik, Introduction to Radar Systems, Second 
Edition (McGraw Hill, New York, 1980). 

[19] A typical SLM would have square pixels of finite di- 
mension and, for computational ghost imaging, be used 
to impress an array of random phases (uniformly dis- 
tributed on [0, 2ir]) on an incident laser pulse. For ana- 
lytical convenience, we are passing to a continuum limit 
for transverse coordinates, rather than discretizing space 
into pixels in these correlation functions. Thus 2ao and 
po roughly correspond to the SLM's diameter and pixel 
size. 

[20] J. H. Shapiro, in Laser Beam Propagation in the At- 
mosphere, edited by J. W. Strohbehn (Springer- Verlag, 
Berlin, 1978), Chap. 6. 

[21] J. W. Goodman, Speckle Phenomena in Optics: Theory 
and Applications (Roberts, Englewood, 2007). 

[22] J. H. Shapiro, B. A. Capron, and R. C. Harney, Appl. 
Opt. 20, 3292 (1981). 

[23] The choice of a Gaussian pupil function is made for an- 
alytical convenience; a unity-transmission circular pupil 
could be used instead at the cost of more complicated 
formulas for noise terms in the SNR analysis to follow, 
as well as in the performance comparison with the laser 
radar. 

[24] G. R. Osche, Optical Detection Theory for Laser Appli- 
cations (Wiley-Interscience, Hoboken, 2002). 

[25] With this assumption, T' v appearing in our spatial- 
resolution formula reduces to T p . 

[26] Because the pulse-repetition interval, T p , is fixed, letting 
N — > oo will necessarily drive the time-averaged turbu- 
lence factor 7 to zero, hence it does not appear in the 
saturation SNR expression. 

[27] R. E. Meyers and K. S. Deacon, Proc. SPIE 7815, 781501 
(2010). 



